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Dirac operator Hw(mo) cannot appear and near-zero modes are strongly suppressed. By measuring 



We propose a lattice action including unphysical Wilson fermions with a negative mass mo of the 
order of the inverse lattice spacing. With this action, the exact zero mode of the hermitian Wilson- 



the spectral density p{Xw)i we h n d a gap near \\y = on the configurations generated with the 
standard and improved gauge actions. This gap provides a necessary condition for the proof of the 
exponential locality of the overlap-Dirac operator by Hernandez, Jansen, and Liischer. Since the 



number of near-zero modes is small, the numerical cost to calculate the matrix sign function of 



Hw(mo) is significantly reduced, and the simulation including dynamical overlap fermions becomes 
feasible. We also introduce a pair of twisted mass pseudo-fermions to cancel the unwanted higher 
mode effects of the Wilson fermions. The gauge coupling renormalization due to the additional 
fields is then minimized. The topological charge measured through the index of the overlap-Dirac 
operator is conserved during continuous evolutions of gauge field variables. 
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I. INTRODUCTION 



In the construction of the lattice chiral fermions, the conventional Wilson-Dirac operator 



.Djy still plays a crucial role. In the domain-wall fermion 



3| the lattice Dirac operator 



is nothing but the Wilson-Dirac operator in four dimensions, but the fermion field has 
interactions also with its fifth dimensional neighbors. The overlap-Dirac operator 4j D 
contains a matrix sign function of the hermitian Wilson-Dirac operator Hw = j^Dyy as 

D = -[l+ l5 sgn(aH w )}, a = (1) 

(X X \~ s 

An important difference from the usual Wilson fermion is that the mass term is given with a 
negative value mo = — (1 + s)/a of the order of the inverse lattice spacing 1/a. In the limit of 
vanishing gauge coupling, the parameter s must be between —1 and 1 in order to obtain the 
single flavor massless fermion; the sign function is well defined since there is a lower bound 
in the eigenvalue spectrum of \H W \. In the presence of the gauge interaction, H w could 
develop zero or near-zero eigenvalues, which makes the the matrix sign function singular. 
The near-zero mode does not exist for sufficiently smooth background gauge fields {U^(x)} 
satisfying a condition ||1 — P, iu {x)\\ < e for all plaquette variables P^x) with e a small 
number less than ~ 1/20.49 [5]. In the actual numerical simulations, however, ||1 — P /iJ/ (a;)|| 
is larger by an order of magnitude, and the near-zero modes appear quite frequently. 

The origin of the near- zero modes is understood as a local lump of the background gauge 
field (or so-called the dislocation). An analytic example of such a gauge configuration and 
its associated exact zero mode is given in Because such a zero mode is localized in space- 
time, the number of the near-zero modes increases as the lattice volume V is increased. In 
other words, the spectral density pLW) of Hw is non-zero at Xw = 0, which is true at 
any finite value of gauge coupling [3]. The localization property of the near- zero modes 
has recently been studied extensively (si, [9], [hJ, and it is found that they are exponentially 
localized unless one enters the Aoki phase, where the flavor-parity symmetry is spontaneously 



broken 11] . Since the radius of the exponential fall-off is of the order of lattice spacing a, 
the effect of the near-zero modes disappears in the continuum limit, and therefore is a lattice 
artifact. 

The effect of the near-zero modes appears as a small residual breaking of chiral symmetry 
in the domain- wall formulation 46] . Namely, the four-dimensional fermion mode receives 
additive mass renormalization m res when the lattice extent in the fifth dimension L 5 is finite 







The problem is not just that m res is finite, but the suppression is only by 1/L 5 rather 



than by exp(— cL^) as expected for the extended non-zero modes 10, ll4|. For the overlap 
fermion, the residual mass can be made arbitrarily small by projecting out the near-zero 
mode and treating them exactly when one calculates the matrix sign function. The problem 
however manifests itself in the locality of the overlap-Dirac operator. The locality is proved 
only when there exists a gap in the spectral density of Hw near zero 15|. Therefore, the 
existence of the near-zero mode persisting in the infinite volume limit could potentially spoil 
the locality of the overlap-Dirac operator and thus the controlled continuum limit of the 
overlap fermion. Even if there exists non-zero density of the near-zero modes, the locality 
may still be maintained [&], provided that the near-zero modes are exponentially localized 
and its exponential fall-off is sufficiently fast. The effect of the near-zero modes becomes 
irrelevant in the continuum limit, as the localization length scales as the lattice spacing a. 
In this case the locality of the near-zero modes should always be checked in order to make 
sure that the Aoki phase is not entered. 

In this paper we propose the use of the lattice action with two flavors of extra Wilson 
fermions with the large negative mass tuq. Since the occurrence of the near-zero modes 
is suppressed by the fermion determinant, we expect that the spectral density p(Xw), as 
defined in ((7j), vanishes at \w = and the near- zero modes are strongly suppressed. By 
suppressing the near-zero modes, the problem of locality of the overlap-Dirac operator is 
essentially solved. Furthermore, the numerical cost for applying the overlap-Dirac operator 
is significantly reduced, because the cost for projecting out the low-lying eigenmodes is 
proportional to the eigenvalue density of Hw- This is especially important when one wants 
to include the dynamical effect of the overlap fermions in the Monte Carlo simulations. 

The idea of addin g th e extra Wilson fermions is very simple and in fact has been around 
for several years 16,|l7, Q, but detailed numerical study has been missing until recently. 
A preliminary report of this work was included in 191 ]. and a paper by Vranas 20[] was 
submitted very recently, after this work had been essentially completed. 

With the action that forbids the occurrence of the zero mode, the global topological 
charge of the gauge field configuration cannot change through continuous deformation of 
the gauge variable. This is because the spectral flow of Hw{fn) can never cross zero at 
m = m under the continuous deformation. Here we use the index of the overlap-Dirac 
operator constructed from Hw{mo) as a definition of the topological charge. Unlike the 
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measurement of the FF operator after some cooling, this always gives an unique integer for 
the topological charge. The property of conserving topological charge resembles that of the 
continuum theory. In the continuum theory, there is an infinitely high barrier in the action 
between different topological sectors, but the barrier is lowered by the lattice regularization 
of gauge field. With the extra Wilson fermions the infinite potential barrier is recovered at 
any finite lattice spacing, and thus the continuum gauge field is better approximated. 

For the Monte Carlo simulations including the fermion determinant, the Hybrid Monte 
Carlo (HMC) 21( provides the most efficient algorithm. HMC is based on the molecular dy- 
namics evolution of gauge field variables under the Hamiltonian including a pseudo-fermion 
action x^H^X- The gauge field evolves with small time steps, that approximate the con- 
tinuous evolution. When a near-zero mode of H\y appears and further approaches the zero 
point, the pseudo-fermion action increases rapidly and generates a repulsive force. This 
means that the topological charge cannot change during the HMC simulation, as far as the 
step size is taken small enough. If the step size is so large that the transmission through 
the potential barrier is allowed, the conservation of the Hamiltonian becomes poor and the 
acceptance rate of the Monte Carlo would become very low. 

Since the topological charge does not change, the HMC simulation is confined in an ini- 
tially given topological sector, and the correct sampling of the 6 = (or any finite value 
of 9) is not possible. This is an advantage for the simulations in the e-regime of the chiral 
perturbation theory (ChPT) 22j, for which one needs the gauge ensemble in a given topo- 
logical sector. For this purpose, some of the present authors tested a modified plaquette 



gauge action to suppress the topology change 



231 ] (see also 



24] ) . The proposal of the present 



work is more solid, as it strictly prohibits the topology change. Out of the e-regime, on the 
other hand, the fixed topology is a disadvantage. But the possible errors due to the incorrect 
sampling of the 6 = vacuum disappear quickly as the physical volume is increased. 

Although we aim at performing dynamical fermion simulations using the overlap formal- 
ism, this paper focuses on the quenched study. Namely, the dynamical overlap fermion is 
switched off, while the extra Wilson fermions are treated dynamically. We then numerically 
study the spectral density of H w for various choices of gauge actions. 

The rest of this paper is organized as follows. In Section [Til we introduce the lattice 
actions that we studied in this work. An implementation in the HMC simulation is also 
presented. Our numerical simulations are explained in Section IIHI Section IIVI contains the 
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main results of this work, i.e. the spectral density of Hw with and without the extra Wilson 
fermions. The topology conservation is a key property of the extra Wilson fermions. We 
discuss how it works (Section |Vj and what is its effect on physical quantities (Section IVip . In 
Section IVII[ we calculate how much the gauge coupling is renormalized by the extra Wilson 
fermions, and show that the finite renormalization can be made small by further adding 
pseudo-fermions with a twisted mass term. Section IVIIII contains our conclusion as well as 
some discussions about possible artifacts due to the fixed topology. 



II. LATTICE ACTION AND ITS IMPLEMENTATION 



We consider a lattice action 

S = S G + S E , (2) 

where Sq is any gauge action, such as the Wilson, Liischer-Weisz, Iwasaki, etc., and Se 
denotes the extra Wilson fermion. In the following numerical analysis, we use for the gauge 
part Sq the plaquette gauge action Spi, with and without a modification to suppress dis- 
locations, and the Renormalization Group (RG) improved (or Iwasaki) action Srq. The 
plaquette action is written as 

f £ £ i n~^T P p?(^/ ' When 1 - ReTri V(*)/3 < e, 
Spi=l x^u 1 - (l-RelVP^(:z)/3)/e ; ( 3 ) 

[ oo otherwise 
where P^fa) is the plaquette variable at x on the \x-v plane. The denominator is introduced 
;o s upp ress the local lump of the gauge field for which 1 — ReTrP /t „(a;)/3 takes a large value 
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25] . When 1/e = 0, it reduces to the standard Wilson gauge action. For this gauge 
action with a finite 1/e, 1 — ReTrP^a;) /3 cannot become larger than the parameter e. If e is 
chosen smaller than 1/20.49, then the hermitian Wilson-Dirac operator is proved to have a 
gap [5J], but we take a larger number 3/2 in the numerical simulations, because otherwise the 
lattice spacing becomes too small even at the strong coupling limit. Another choice of the 
gauge action is that of Iwasaki [26] , which includes the rectangular term with the parameter 
C\ = —0.331. The denominator as in ([3]) is not introduced for this case, but it is known that 
the rectangular term has an effect to suppress the dislocations. 
The extra fermion term Se is written as 

s e = ^2jf>(x)Dw(mo)if3(x) +^2(j)\x)[D w (m ) +^7 5 r 3 ]0(x), (4) 



where ip denotes two flavors of extra heavy Wilson fermion with a negative mass rriQ. The 
second term is a pseudo-fermion term introduced in order to cancel unwanted effects of the 
Wilson fermion especially in the ultraviolet region, which leads to a large shift of the (3 value 
to be used in the simulation as discussed later in Section IVHI Because of an additional mass 
term (pHu'y^r^, which is twisted in the flavor space by r 3 , the extra Wilson fermion works 
to suppress the near-zero modes of Hw(rn ) as expected. In fact, the action Se generates 
the suppression factor 



det ""T;» 5 (5) 
_H w (m y + [i\ 

in the partition function. The twisted mass fi controls the range of the near-zero eigenvalues 
suppressed by the numerator. The eigenvalues whose absolute value is lower than \i are 
strongly suppressed, while the other higher modes are less affected. The limit of \i — > 00 
corresponds to the case where the pseudo-fermion term is switched off. When /i = 0, the 
cancellation is exact, and there is no extra fermions and pseudo-fermions. 

Since the action (j2J) includes the fermion, some dynamical fermion algorithm is needed 
to generate the gauge field ensembles. Application of the HMC algorithm is straightforward 
except for the additional boson term. In order to cancel the higher modes of H\y efficiently, 
we use only one pseudo-fermion for both fields. Namely, the Hamiltonian for the molecular 
dynamics evolution contains a term 

E^Wfl^H + ^75][^ ) H/("lo)] _1 [^( m o)] _1 [^ , H/( m o) "WlJxW (6) 

x 

with the pseudo-fermion field \. Then, the fermion force derived from (jfJJ) largely cancels in 
the combination [Dw(m ) + 2^75] [Dw(mo)] -1 , when the twisted mass \i is small. 

In the molecular dynamics evolution with ([6]), an inversion of the Wilson-fermion matrix 
D\ v {mo)Dw{fT r Lo) is necessary, and therefore it costs much more than the usual gauge action. 
The cost is, however, not substantial compared with the inversion of the overlap-Dirac 
operator, for instance, needed for the dynamical overlap fermion simulation. The cost for 
the inversion of is proportional to the inverse of lowest-lying eigenvalue A^ in , which 
is lifted by the introduction of the suppression factor (J5j). It means that the cost does 
not increase arbitrarily, even though there is no explicit lower limit on the lowest-lying 
eigenvalue. 



6 



III. HMC SIMULATIONS 



We performed Monte Carlo simulations including the extra Wilson fermions. Although 
the fermions are included, they are unphysical and irrelevant in the continuum limit. The 
physical sea quarks are not included. 

The numerical simulations have been done on a 16 3 x 32 lattice with three choices of the 
gauge actions. 

1. Spi with 1/e = 0, the standard Wilson gauge action. 

2. Spi with 1/e = 2/3. With this choice the plaquette variable P^ u (x) can take any 
value in the SU(3) gauge group except for the two points e ±l27T ^ 3 I with / the 3x3 unit 
matrix. At these isolated two points, ReTrP M ^(a;)/3 becomes minimum. The positivity 



violation as argued in [27| occurs for e smaller than 3/2. 

3. Srg, the Iwasaki gauge action. 

The simulation parameters are listed in Table [fl For each of the three choices of the gauge 
action we take three values of //, the twisted mass of the extra pseudo-fermions, to be 0, 0.2, 
and 0.4. The large negative mass mo is always set to —1.6. The gauge coupling j3 is chosen 
such that the lattice spacing determined through the Sommer scale r$ (~ 0.49 fm from 
phenomenological models) is roughly tuned to 0.125 fm, which corresponds to r /a = 3.9. 

All the simulations are carried out using the HMC algorithm, including those with [i = 0, 
i.e. no extra fermions and pseudo-fermions. The step size St of the molecular dynamics 
evolution with the leapfrog scheme is 0.01 for all lattices except for one choice of the action 
(Spi with 1/e = 2/3 and n = 0.4). The unit length of the HMC trajectory is set to 0.5. 
The total number of the HMC trajectories and observed acceptance rate P acc are listed in 
Table [0 Also listed is the integrated autocorrelation time T int measured for the plaquette 
variable. 

It is known that the HMC simulation becomes unstable due to (near-)zero eigenmodes 



(see, for example, 28]). Since the fermion force F in molecular dynamics steps contains 
a piece that is proportional to 1/A m » n , the step size 5t must be kept small such that the 
combination F x 5r is less than 0(1). Because there is no explicit lower bound on X m in 
for the Wilson-type fermions, very small eigenvalue could appear, though it is suppressed 
by the fermion determinant. The problem appears as exceptionally large values of AH, the 
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TABLE I: Simulation parameters. 




traj. traj. 

FIG. 1: History of AH for the lattice "PI" with 1/e = 0. The left panel shows our main run with 
5t = 0.01, while the right panel is for a reduced step size: 5t = 0.005. 

difference of the Hamiltonian between the initial and final steps of the HMC trajectory. In 
Figure [T] we show a typical HMC history of AH for the lattice "PI" with 1/e = (standard 
Wilson gauge action). On the left panel we show our main run listed in Table [H AH is of 
order one for most trajectories but develops several large spikes during long runs. We also 
made a test run with a reduced step size by a factor of 2 (St = 0.005) and show the result on 
the right panel, for which the fluctuation of AH is reduced and the occurence of the spikes 
becomes rare. 

The volume dependence of the computational cost for the extra Wilson fermion is quite 



S 



normal. We repeated the test run at St = 0.005 on a smaller lattice 12 3 x 24, and measured 
the average value of AH. With 1,000 HMC trajectories we obtain (AH) = 0.0038(23), which 
may be compared with the result on a larger lattice 0.0061(59). Although the statistical 
error is large, the volume scaling is consistent with the conventional 0(V) 2t| (or milder), 
which implies that the total simulation cost scales as V 5 / 4 . There is no significant spike on 
both lattice at this small step size. 

The exceptional trajectories result in a violation of area preserving property of the 
leapfrog integration of the molecular dynamics evolution. If it is violated, the detailed 
balance condition of the Monte Carlo algorithm cannot be proved, and therefore the exact- 
ness of the algorithm is lost. This is a problem associated with a discontinuity in the effective 
pseudo-fermion action, as in the case for ([6]). If the step size is small enough, the trajectory 
would never pass the discontinuity because of the strong repulsive force from the potential 
wall, and no problem arises. But for some large step size the trajectory may occasionally go 
across the discontinuity, and then the conservation of the Hamiltonian is strongly violated. 

This problem can be monitored by the quantity (e~ AH ), which must be consistent with 
1 when the area preserving property is satisfied. In the second last column of Table [I] we list 
the value of (e~ AH ). For most cases it is consistent with 1 within one standard deviation. 
The problem is expected to disappear if one chooses small enough step size Sr. In fact, for 
the run "PI" with 1/e = 2/3 and \i = 0.4, for which 5r is 2/3 of other cases, the number of 
large spikes is much reduced and the relation (e~ AH ) = 1 is satisfied to a good precision. 

In the following analysis we assume that the gauge configurations are properly sampled. 
Even if the exceptional trajectories exist, they are almost always rejected by the Metropolis 
test and thus do not affect the following trajectories. In the productive run with the dy- 
namical overlap fermion, that we are currently carrying out, the step size is carefully chosen 
to avoid the potential problem. 



IV. SPECTRAL DENSITY 



We investigate the effect of the extra Wilson fermions on the spectral density of the 
hermitian Wilson-Dirac operator ify^(mo) in the quenched approximation. 

We take 30-100 gauge configurations for each runs listed in Table |T] and calculate the 
low-lying eigenvalues of Hw(mo) with m Q = —1.6. For each runs the gauge coupling is 
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chosen such that the lattice spacing is roughly tuned to 0.125 fm. Statistical correlation 
between consecutive gauge configurations is negligible, as they are separated by 200 HMC 
trajectories. We use the conventional Lanczos algorithm, and calculate the eigenvalues with 
its absolute value less than 0.2 in the lattice unit. 

In Figure [2] we plot the absolute value of the observed eigenvalues \\y for each set of 
configurations. Different choices of lattice actions are shown in separate panels, and in each 
panel the results at three values of p (0.0, 0.2, and 0.4) are plotted. As one can clearly see, 
there are significant number of eigenvalues less than 0.01 for any of the three gauge 
actions if the extra Wilson fermions are not included (p = 0). On the other hand, with the 
finite values of p, the eigenvalues less than 0.02 do not appear at all. The difference between 
^=0.2 and 0.4 is not significant, but p=0.2 seems slightly better. 

These observations can be made more quantitative using the spectral density p(Xw)- The 
spectral density is defined by 



where (• • •) denotes the ensemble average and A n represents each eigenvalue on a given gauge 
configuration. In Figure [3] we plot the histogram of p(Xw)- Each bin in the plots has a size 
of 0.0192 in \w- We observe non-zero spectral density at \w = when there is no extra 
fermion introduced (p = 0.0). We calculate the value of p(0) by fitting an integrated density 
/(A) = f^d\p{\) with a polynomial 2p(0)A + O(A 3 ). The results are p(0) = 1.509(3) x 10~ 3 
for "PI" with 1/e = 0, 0.494(1) x 10~ 3 for "PI" with 1/e = 2/3, and 0.177(1) x 10~ 3 for 
"RG". As well known, the density of the near-zero modes is substantially reduced for the 
improved actions. 

With the extra Wilson fermions the near-zero mode density is essentially zero; there is 
no events in the lowest bin. The curve of p(Xw) near -W = is consistent with ~ A^/ for 
both p = 0.2 and 0.4. The difference between p = 0.2 and 0.4 is not significant, except for 
the highest bin at \Xw\ — 0.2. 

Among the three choices of the gauge actions, the suppression of the near-zero modes 
is most efficient for the RG action. We should note that the scale of the vertical axis in 
Figure [3] is different among the three plots, and the value at |Aw| — 0.2 is 6-10 times smaller 
for the RG action than the standard Wilson gauge action (1/e = 0). The choice "PI" with 
1/e = 2/3 is in between the other two. Therefore, we may conclude that the rectangular 
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FIG. 2: Low- lying eigenvalues of Hw(mo). Lattice actions are "PI" with 1/e = (top), "PI" with 
1/e = 2/3 (middle), and "RG" (bottom). Data for three values of jj, (fi = 0.0, 0.2, and 0.4) are 
shown in each plot. The highest mode is also shown, whose value is ~ 5.9, almost independent of 
the gauge configuration. 
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FIG. 3: Histogram of the spectral density of Hw(mo). Lattice actions are "PI" with 1/e = 
(top), "PI" with 1/e = 2/3 (middle), and "RG" (bottom). Data for three values of y. (p, = 0.0, 0.2, 
and 0.4) are shown in each plot. 

term is more effective to suppress the dislocation than restricting the range of the plaquette 
variable by the denominator in (J3j). 
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V. TOPOLOGY CONSERVATION 



As we discussed in the Introduction, the net topological charge defined as an index of 
the overlap-Dirac operator cannot change by the continuous deformation of the gauge field 
variables when the extra Wilson fermion is included. This is approximately the case in our 
simulation using the HMC algorithm with small molecular dynamics time step St. The 
large spikes in AH as discussed in Section HTT1 may indicate some attempts of the gauge field 
to go beyond the potential barrier, but the acceptance probability of such trajectories is 
essentially zero and therefore the topology does not change in the accepted trajectories. 

In order to explicitly check the topological charge of the generated gauge configurations, 
we calculate the near-zero eigenmodes of the overlap-Dirac operator D as defined in ([T]). 
The kernel operator is the Wilson-Dirac operator with tuq = —1.6, the same value as in 
the extra Wilson fermion. The topological charge can be identified as the number of the 



30] 



left-handed or right-handed zero modes. In practice we adopt the method proposed in 
Namely, we calculate the low-lying eigenvalues of chirally projected operators = P±DP±. 
Since the non-zero eigenvalues appear in pair between D + and D~, the remaining unpaired 
eigenvalues very close to zero can be identified as chiral zero modes. 

We calculate the eigenvalues of D on 80 configurations of "RG" with /i = 0.2, and find no 
exact chiral zero mode. This is consistent with our expectation, as these configurations are 
generated starting from an initial configuration in the trivial topological sector, i.e. Q = 0. 
In addition, we also generated another set of configurations starting from a Q = —2 config- 
uration (number of samples is 60). On all of these configurations we confirm the presence 
of two left-handed zero modes as expected. From these observations we conclude that the 
topological charge is indeed conserved during the HMC simulations. We are currently per- 
forming many set of simulations including dynamical overlap fermions as well as the extra 
Wilson fermions. Also in these simulations, there has been no sign of the change of the 
topological charge so far. 

VI. TOPOLOGY DEPENDENCE 

Since the topological charge conserves during the HMC simulation, one cannot sample 
the correct 9 = vacuum, for which the topological charge distributes according to the 
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topological susceptibility. However, the cluster decomposition property of the local field 
theory suggests that the physical quantities measured in a local space-time region do not 
depend on the topological fluctuations occuring far apart from that region. This means that 
the physical quantities do not depend on the global topological charge, as far as the space- 
time volume is large enough. The topological fluctuation is controlled by the topological 
susceptibility Xt = (Q 2 )/V, which is proportional to mS = M 2 F 2 in unquenched QCD near 
the chiral limit. The typical length scale for the topological fluctuation is then given by 
(M n F n )-^ 2 . The lightest possible pion mass in our planned dynamical simulation is M w ~ 
300 MeV, for which the typical scale is ~ 1 fm. It becomes even smaller for heavier sea 
quark masses or for the quenched theory. Therefore, our naive expectation is that the fixed- 
topology simulations in a (2 fm) 4 box or larger do not have too much impact on physical 
quantities. 

On a rather general ground, one can prove that the hadron masses in a fixed topological 
sector deviate from the correct value in the 6 = vacuum by 0(1/ (Q 2 )), which scales as 
~ 1/V 3J|. An estimate of the coefficient using the chiral effective lagrangian implies that 
the deviation of pion mass does not exceed 1% even for light pions (~ 300 MeV) on a (2 fm) 4 
box, and the effects on other hadrons are even smaller. 

In this paper we show a study in the quenched approximation, instead of the time- 
consuming dynamical fermion simulations. In the quenched theory, however, the situation 
could be drastically different because the hadron masses are more infra-red sensitive due to 
the double pion pole structure at the one-loop level of quenched chiral perturbation theory. 
Because of the integral of the form Jd 4 pml/(p 2 + m 2 ) 2 , the diagram is infra-red divergent 
in the chiral limit, and the information of the topological fluctuation of the whole space- 
time is gathered through the singlet vertex tUq proportional to the topological susceptibility 
Xt, whose integral over space-time gives Q 2 the fixed topological charge in this case. We 
therefore expect the quenched chiral logarithm that depends on the topological charge as 
Q 2 /V rather than its average (Q 2 )/V. 

In order to see the Q dependence in the simulations, we carry out a calculation of pion 
mass on quenched gauge configurations with fixed topological charges \Q\ = 0, 2, 3, and 
6. The gauge configurations are generated with the extra Wilson fermions at (3 = 2.37 
(RG gauge action), and the statistics is 100 for each topological charge. The results for 
(amps) 2 /(am) are plotted in Figure H] as a function of the quark mass am. We find rather 
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FIG. 4: Pseudo-scalar mass squared calculated on the quenched configurations at fixed topological 
charges \Q\ = 0, 2, 3, and 6. 

large topological charge dependence as expected in the quenched theory: the difference of 
amps between \Q\ = and 6 is as large as 10%. A similar dependence is also found in [32 ]. 

Such a large topology dependence is not expected for unquenched QCD, because there 
is no pathological infrared divergence. In our on- going dynamical overlap fermion project, 
we are planning to investigate the topology dependence in detail as well as the volume 
dependence. 



VII. BETA SHIFT 



The extra Wilson fermions and their associated pseudo-fermions have masses of order of 
the lattice cutoff and do not affect the physics at low energy. Their only effects are ultraviolet 
ones, such as the renormalization of the strong coupling constant and quark masses. Here 
we calculate the finite renormalization of the strong coupling constant (or the j3 shift) due 
to the extra fields at the one-loop level in perturbation theory. Such calculation for the 
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massless Wilson fermion was done before 
the large negative mass. 

The shift of the coupling constant is given as 

1 



341 ]. but we must repeat the calculation with 
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where II (p) is the vacuum polarization function of gluon due to fermions: 
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FIG. 5: Tl{p)/Nf for various mass parameters. 



For the usual (lattice) fermions, it is written as 

1 



n(p) = n 



247T 2 



ln(aV) + 



(10) 



and the constant kt depends on the fermion formulations. One- loop calculation gives 



—0.013732 (Wilson fermion) 35j, |36|, —0.038529 (clover fermion) [37]. The logarithmic 



term in (jTOl) corresponds to the usual one-loop running of the coupling constant and gives 
the Nf dependent coefficient of the beta function b = (jjN — |-A^/j. 

The unphysical fermions with a large negative mass do not contribute to the logarithmic 
term in (jTOl) . Therefore, the p — > limit is finite. Figure [5] shows H(p)/Nf for the negative 
mass Wilson fermion with masses am = —0.80 ~ —1.80. Since the numerical integral 
becomes unstable near ap = 0, we obtain the ap — > limit by an extrapolation using the 
data points below ap = 1. The numerical results for 11(0) are listed in Table HT1 

At am = —1.6 the (3 shift 5(3 = —611(0) is calculated as 0.786(4) at the one-loop level. 
It means that we have to reduce the (3 value by this amount in order to simulate without 
changing the lattice spacing. The actual value could be significantly larger due to higher 
order corrections. With the mean field improvement, for instance, the value is divided by 
the expectation value of the plaquette (P), which is typically around 0.6 for the Wilson 
gauge action at (3 = 6. 

Including the twisted mass pseudo-fermions, the (3 shift is substantially reduced. It 
vanishes in the limit of \i — > and the contribution for finite \x starts from fi 2 . Results are 
shown in Figure [6]for am = —1.6 and the twisted mass \x = 0.1-0.4. Note that the vertical 
axis is one order of magnitude smaller than that in Figure [5j At p, = 0.2, the f3 shift is 
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FIG. 6: H(p)/Nf including the twisted mass pseudo-fermions. 
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-0.00844(7) 


-0.0133(1) 


-1.8 


-0.0952(4) 




-0.01209(4) 







TABLE II: Vacuum polarization function U(p)/Nf at p = 0. 



0.0502(4) (N f = 2), which may be compared with the actual data: f3 = 5.83 -> 5.70 for "PI" 
with l/e=0, 2.33 -> 2.23 for "PI" with l/e=2/3, and 2.43 -> 2.37 for "RG" . The one-loop 
calculation underestimates the measured value by about a factor of two for the plaquette 
gauge action, while it gives a good approximation for the RG action. This is consistent 
with the expected mean field enhancement of the coupling constant for the plaquette gauge 
action. 

The small (3 shift for the cases with the twisted mass pseudo-fermions is desirable, because 
one can avoid too small (3 values in the dynamical fermion simulations. In particular, for the 
Wilson gauge action there is a remnant of the fundamental-adjoint phase transition 38] in 
the strong coupling regime ((3 < 5.2). With dynamical fermions it may appear as a real first- 
order phase transition, which prevents one from taking smooth continuum limit. Examples 
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are found in the simulations with two 39|, |40| and three flavors 4l| of the Wilson-type 
fermions. 



VIII. CONCLUDING REMARKS 

By introducing the extra Wilson fermions with the large negative mass, it is possible to 
provide a gap in the spectral density of H w and thus to remove the problem associated 
with the near-zero modes. For the overlap fermion the exponential localization of the Dirac 
operator is guaranteed; for the domain-wall fermion the residual chiral symmetry breaking 
will be much reduced. For three choices of gauge actions, we confirm that the spectral density 
of Hw vanishes at Xw = and the occurrence of the near- zero mode is very much suppressed. 
The number of remaining near-zero modes is smallest for the Iwasaki gauge action. By also 
adding the twisted mass pseudo-fermions, the effect on the coupling renormalization can be 
minimized while keeping the good property of suppressing the near-zero modes. 

A remaining problem concerning the locality of the overlap-Dirac operator is its actual 
localization length in the numerical simulations. For the numerical simulation with con- 
trolled discretization error, the localization length has to be much smaller than 1/Aqcd, 
and this gives a stronger practical upper limit of the gauge coupling in addition to the more 
fundamental limit from locality. Golterman, Shamir, and Svetitsky argued that the local- 
ization length is controled by the mobility edge of the eigenmodes of Hw |s[ O, Q. We are 
currently calculating the spatial correlation on the low-lying mode eigenvectors in order to 
determine the mobility edge. This work will be published elsewhere. 

Although the extra Wilson fermion is also useful in the quenched QCD simulations, our 
main objective is to perform the unquenched simulations using the fermion formulation with 



exact chiral symmetry. So 



only on small lattices 42 



ar, the dynamical overlap fermion simulation has been attempted 



43 



44 



451 ]. as its computational cost is extremely large. One of 
the reasons for the large numerical cost is the treatment of the sign function. When the low- 
lying eigenvalue Hw passes zero during the molecular dynamics evolutions, the sign function 



changes its value discontinuously and its derivative diverges. Fodor, Katz and Szabo [42] 
introduced so-called the reflection/refraction process on the Xw = surface, which requires 
monitoring of the low-lying eigenvalues and takes lots of computational costs. The cost 
would scale as V 2 , as the lattice volume V is increased, and therefore the simulation on 
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large lattices could be prohibitively costly. This problem can be totally avoided by the 
extra Wilson fermion, because the zero crossing never happens if the molecular dynamics 
step size is chosen small enough. The extra cost for the Wilson fermion is negligible in the 
dynamical overlap fermion simulation. Using this lattice action, we are currently carrying 
out the dynamical overlap fermion simulations on a 16 3 x 32 lattice. 

An immediate question on the strategy of suppressing the low-lying modes of Hw is 
that the simulation is trapped in a given topological sector and one cannot achieve the 
correct sampling of the 9 = vacuum of QCD. This is an algorithmic problem of the 
molecular dynamics simulation, in which the gauge variables are changed continuously along 
trajectories. We first emphasize that the continuum theory has the same property, i.e. the 
topological nature of the gauge field on a torus. Any lattice gauge action should recover 
this topological property as the continuum limit is approached, and our choice has this 
property at any finite lattice spacing. The configuration space of a given fixed topology is 
simply connected in the continuum limit, and therefore the ergodicity of the Monte Carlo 
simulation is satisfied even with a fixed topology. At finite lattice spacings there is no proof 
of the ergodicity, but there is no indication of its violation either. In order to detect such 
an effect if any, we should probably look at quantities sensitive to the topological charge 
density. That is an interesting subject, which we leave for future studies. 

Second, because of the cluster decomposition principle any physical quantity should not 
depend on the net topological charge, if the lattice volume is large enough. Even on a 
lattice with a fixed topological charge, local topological fluctuation could occur at different 
positions on the lattice keeping the net topological charge constant, and the frequency of 
the topological fluctuation is characterized by the topological susceptibility. In other words, 
the effect of fixing the net topology is a finite volume effect. (An important exception to 
this statement is the CP-odd quantities, such as the neutron electric dipole moment.) As 
we discussed in Section [VTl hadron masses calculated at a fixed topology contains a finite 
volume correction of 0(1/ V) to the physical value in the 9 = vacuum 31]. Although the 
effect is quite large in the quenched QCD due to the sickness of quenched theory, we expect 
a small effect on dynamical lattices. Study of the topological charge dependence is underway 
as a part of our project of dynamical overlap fermion simulations. 
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